Anomalous Scaling of the SO (8) Symmetric Phases in the Two-Leg Ladder 
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We carry out a complete analytical stability study for the SO (8) symmetric phases in the weakly- 
interacting two-leg ladder. It is shown that the SO(8) symmetry is robust under generic perturba- 
tions. Since there are no fixed points in the one-loop renormalization group equations, the conven- 
tional classification of relevant and irrelevant perturbations fails in this case. A new classification 
is defined and explained in detail. It leads to the anomalous scaling in the ratios of excitation gaps 
and an universal exponent 1/3 is extracted. This new method is also applied to the well studied 2D 
Ising model and similar exponent is calculated. 



Introduction. — Ladder materials have attracted lots of 
attentions in the past decade because of both theoretical 
and experimental interest Q. In particular, lots of efforts 
were focused on the simplest two-leg ladder [p|-|il"f . With- 
out doping, there is one electron per site on average (usu- 
ally referred as "half-filled" ) . Because of the Coulomb 
repulsion, charge excitations acquire a gap which makes 
the two-leg ladder a Mott's insulator. Besides, due to the 
tendency for singlet bond formation across the rungs of 
the ladder, spin-liquid behavior is expected. Both numer- 
ical and analytical approaches support that the ground 
state at half filling is in the Mott-insulating spin-liquid 
phase @-g]. 

In general, one expects that charge and spin gaps could 
be rather different. Surprisingly, a complete degeneracy 
of charge, spin and single particle gaps emerges in the 
asymptotic weak coupling limit. Based on the one- loop 
renormalization group (RG) analysis, the effective low- 
energy theory for the two-leg ladder is described by the 
exactly soluble Gross-Neveu model with a global SO (8) 
symmetry. Other than the exact degeneracy of excitation 
gaps, many interesting results can be drawn from the 
Bethe Ansatz solution. There have been critics [ ^2|Jl3| 
about the unexpected restoration of the SO (8) symmetry 
in weak coupling, concerning about the stability of the 
symmetry under generic perturbations and the limitation 
of the RG equations derived in weak coupling. It was 
argued that the SO (8) symmetry derived by perturbative 
calculations is fragile because the system eventually flows 
toward the fixed point in strong coupling. 

In this Letter, we address both points by performing 
a complete an analytical study for the stability in the 
vicinity of the SO (8) symmetric phases. It should be 
emphasized that this beautiful symmetry restoration in 
weak coupling does not implies that the fixed point in 
strong coupling is SO (8) symmetric at all. On the con- 
trary, the symmetry reflects the local topology of RG 
flows near the trivial Fermi liquid fixed point which can 
be safely described by perturbative calculations. Indeed 
it is rather obvious that, when the interaction U is much 
larger than the bandwidth f, the charge gap is much 
larger than the spin gap. Thus, the symmetry in the 
strong coupling is completely destroyed. As to the sta- 



bility of the symmetry under generic perturbations, the 
complete stability check shows that the SO (8) symmetry 
is robust in weak coupling. Even though some couplings 
are relevant according to the conventional classification, 
they do not destroy the symmetry but only give rise to 
anomalous corrections which scale as (CZ/i) 1 / 3 . A new 
classification of "relevant" and "irrelevant" perturbations 
is necessary here. Finally, we apply the new method to 
the well-known 2D Ising model and show that the anoma- 
lous scaling behavior near the critical point can be very 
general, as long as more than two relevant couplings are 
present. 

The SO (8) symmetric rays. — For the two-leg ladder at 
half filling, the number of possible interactions is greatly 
reduced to nine in weak coupling. Within the one-loop 
RG calculations, these nine couplings gi are described by 
a set of coupled first-order differential equations, 



dl 



(1) 



where A^ k are nine 9x9 constant matrices. These ma- 
trices can be found by rewriting the RG equations given 
in Ref. (9j. For a generic interacting Hamiltonian, the 
bare values for these nine couplings can be straightfor- 
wardly determined. However, it is generally very difficult 
to obtain the solution for these coupled flow equations in 
analytical form. 

Simple analytical solutions emerge if the interactions 
are chosen in a specific way. These special solutions 
are later referred as "symmetric rays" . Suppose the 
bare couplings of the specific interacting Hamiltonian are 
gi{0) = i"ig(0), where g(Q) — (U/t) <C 1 is small while r$ 
are order one constants which satisfy the algebraic con- 
straint. 



n = A 



jk r j r k- 



(2) 



It is straightforward to show that the ratios between 
couplings remain the same and the nine complicated 
equations reduce to single one, g = g 2 ! The solution 
is g(l) = l/(ld — I), where the divergent length scale 
U = (t/U). Of course, one should keep in mind that 
the solution g(l) is only valid when it does not flow out 
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of the weak coupling regime. These special Hamiltoni- 
ans, whose couplings are described by these symmetric 
rays, turn out to be SO (8) symmetric. For the two-leg 
ladder, four different phases (named as D-Mott, S-Mott, 
CDW and SP in Rcf. |]) are of the central concerns. 

It was shown previously that the two-leg ladder in 
weak coupling always scales into one of the four different 
symmetric phases 0. However, this numerical approach 
was criticized that the SO (8) symmetric phases might 
have instabilities which happen not to be tackled by the 
limited types of interactions considered in the numerical 
study. To make up the fissure, a complete stability check 
near the SO (8) symmetric rays is desirable. 

Stability analysis. — To describe the RG flows in the 
vicinity of the SO (8) symmetric rays, it is sufficient to 
consider the linearized version of Eq. ([!]) . For a generic 
interaction, the couplings are separated into symmetric 
and asymmetric parts, gi(l) — fig{l) + Agi{l). In the 
vicinity of the symmetric rays, the deviations are small, 
A<7i (Z) <C g{l). Keeping the leading order term, the lin- 
earized RG equations are 



d(Ag t 
dl 



B, 



{Id - I) 



A 3l 



(3) 



where B{j — 2A l - k rk. The matrix By can be brought 
into diagonal form by a linear transformation. As a con- 
sequence, the RG equations decouple into nine indepen- 
dent ones, 



dl 



A; 



(Id - I) 



5gi 



(4) 



where 5gi are couplings after the linear transformation 
and Xi are the eigenvalues of the matrix Bij. Although 
the matrix Bij are different for each SO (8) symmetric 
rays, the eigenvalues are identically the same 



Ai = 2, 



2 2 2 



3' 3' 3' 



(5) 



This coincidence implies that the results of the stability 
check only rely on the symmetry group but not on the 
details of the phases [ jl4j . 

So far, we have a single equation, g = g 2 , describing 
the renormalization along the symmetric rays and nine 
for deviations from the rays as in Eq. (^). Apparently, 
there must be one redundant equation among them be- 
cause the original number of equations is only nine. It 
doesn't take long to find out that it corresponds to the 
flow equation with largest eigenvalue A = 2. This cor- 
responds to a trivial case that all couplings are shifted 
along the symmetric ray, i.e. 8gi — ri8g. Since the flow 
along the ray is described by g — g 2 , linearization leads to 
5g = 2g{l)Sg. Since we only need eight equations to de- 
scribe the deviations from the symmetric rays, the A = 2 
equation is only an artifact and should be ignored. 



The other eight eigenvalues describe how the flow goes 
once the bare couplings are off the ray. Starting from 
the bare values 5<?j(0) = 5rig(Q), where 5i\ <c r,*. The 
solutions of Eq. (0) are 



5 9i (l) 



{i d - iy 



l-A, 



(6) 



According to conventional classification jlj], for A.; < 
the deviations diminish under RG transformations as 
shown in FIG. 1(a) and thus are classified as irrelevant 
couplings. For Ai > 0, any small deviations get enhanced 
and the flow is pushed away from the symmetric rays as 
in FIG. 1(b) and 1(c). These are classified as relevant 
couplings. The conventional wisdom tells us that there 
are three relevant couplings (Aj = 2/3) and five irrelevant 
ones (Ai = —1/3). Since a generic interaction in princi- 
ple could generate asymmetric deviations in all couplings, 
one might rush to the incorrect conclusion that the SO (8) 
symmetry is not stable. However, the first-glance guess 
is wrong because the conventional classification is based 
on the perturbative analysis near a fixedpoint while we 
are dealing with running symmetric rays. A new set of 
rules to identify relevant perturbations is in order. 



(a)X<0 (b)0<X<l {c)X>l 

FIG. 1. The topology of RG flows near the symmetric ray 
with (a) A < 0, (b) < A < 1 and (c) A > 1. It is clear that 
the coupling is irrelevant for A < and relevant for A > 1. 
The RG flow is more subtle for < A < 1. In this case, 
although the deviation from the symmetric ray is growing, 
the slope remains the same as the ray. 

The crucial criterion is whether the deviations Sgi{l) 
grow larger than the symmetric coupling g{l) = {1^—1) . 
Just growing larger than the bare values under RG trans- 
formations is not qualified as a relevant perturbation. 
This subtle but important difference is best illustrated by 
calculating gap functions using scaling arguments. If the 
degeneracy of excitation gaps is maintained, the SO (8) 
symmetry is robust and vice versa. 

Anomalous scaling. — Since the RG equations are only 
valid in weak coupling, we cut off the RG procedure when 
the symmetric coupling g{l c ) = 1- At this cutoff length 
scale l c , the deviations in Eq. (n) are 



l-A; 



(7) 



For A.; > 1, the deviations are larger than the symmetric 
coupling and should be classified as "relevant". As long 
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as Xi < 1, the deviations at the cutoff length scale are still 
vanishingly small and should be viewed as "irrelevant". 
This new classification is different from the conventional 
one for < A; < 1. We would see clearly soon that the 
new classification is appropriate for stability check near 
a running symmetric ray. The grey regime < Xi < 
1 between the new and conventional rules gives rise to 
anomalous scaling exponents. 

Under RG transformations, the gap functions scale 
like, 



A, [ff(0) 1 ^(0)] = e- , "A i [l,S gi (l c )], 



(8) 



where Sgi(l c ) are at most of order (U/t) 1 / 3 . The key 
point is that, although some perturbations are enhanced 
to order (U/t) 1 / 3 , which is larger than the bare order U/t 
values, they are still small at the cutoff length scale. The 
effective Hamiltonian can be separated into two parts 
H = Ho + SH - the SO (8) symmetric and asymmetric 
parts. Since the asymmetric part is of order (U/t) 1 / 3 , 
the changes of the gaps would be of the same order by 
standard perturbation theory. Without the deviations, 
the SO (8) symmetry guarantees the exact degeneracy of 
all gaps, i.e. A[l, 0] = A. The presence of perturbations 
modifies the gap functions, 



A* [l,fo(ic)] 



f 



7 1 + 



(9) 



It is clear that, in the weak coupling limit U/t — > 0, the 
degeneracy of all gaps is recovered. It implies that the 
SO (8) symmetry is indeed robust under generic pertur- 
bations. 
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FIG. 2. Ratios of charge, spin and single particle gaps plot- 
ted versus the interaction strength U /t. The ratios approach 
unity in the asymptotic U/t — > limit with anomalous cor- 
rections of order (U/t) 1 / 3 . 

The anomalous scaling exponent 1/3 is clearly seen in 
the ratios between charge, spin and single particle gaps 
(see FIG. 2), 



Ai 

-T— « 1 

A., 



U 

'7 



(10) 



This simple exponent is rather fascinating because the 
RG equations we started from are very messy. But, de- 
spite of which SO (8) symmetric phases the system flows 



into, the exponent of the gap function corrections is uni- 
versally equal to 1 — Xi = 1/3! 




FIG. 3. The RG flows near the critical point in the 2D 
Ising model. Notice that the local topology of flows near the 
critical point is similar to FIG. 1(b). 

2D Ising model. — The anomalous scaling discussed 
here is not limited to the particular two-leg ladder at 
all. As long as more than two relevant couplings (ac- 
cording to the conventional classification) are present, 
this interesting phenomena shows up somewhere. Here 
we use the 2D Ising model as another example. The RG 
flow diagram for the model is shown in FIG. 3. There 
are two relevant perturbations near the critical point - 
the magnetic field h and the temperature deviation from 
the critical point t. It is well known that the magneti- 
zation m[h,t] scales differently with respect to these two 
relevant perturbations, 



m[u, 0] 
to[0, u] 



(11) 
(12) 



Here 1/5 = 1/15 and j3 — 1/8. Although both are rel- 
evant, the magnetic field grow faster than the later un- 
der RG transformations. We would find out that the 
scaling of magnetization m[h,t] in the presence of both 
perturbations looks similar to the scaling behavior with 
only non-zero magnetic field. Suppose now we start 
from a bare coupling in the vicinity of the critical point, 
(h,t) = u(cos9,sin9) with < 6 < n/2. The scaling 
function of the magnetization is m[h,t] — h 1 / 5 F(t A /h) 
with the exponent A = /3<5 = 15/8>1. Since the argu- 
ment inside the scaling function t A /h ~ u A_1 is small, 
we can expand the function around zero. After a bit 
algebra, the magnetization is 



m[u cos 9, u sin 9] 
m [u, 0] 



l+c(9)u 



7/8 



(13) 



Notice that the exponent comes from A — 1 = 7/8. It 
is clear that, when close to the critical point u — ► 0, the 
magnetization scales as if the temperature deviation is 
not present at all. In this sense, the temperature devi- 
ation t should be called "irrelevant". However, it does 
lead to a non-trivial correction with anomalous exponent 
7/8. 

In summary, we have shown analytically that the 
SO (8) symmetric phases are stable in weak coupling. 
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Since we arc dealing with symmetric rays but not fixed 
points, a new classification is denned and explained in 
detail. Corrections to the SO (8) symmetry at finite but 
still weak couplings acquire a non-trivial exponent 1 /3 in 
the ratios of excitation gaps. 
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